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Abstract 

Thermalization of classical fields is investigated in a 4 scalar field 
theory in 1 + 1 dimensions, discretized on a lattice. We numeri- 
cally integrate the classical equations of motion using initial condi- 
tions sampled from various nonequilibrium probability distributions. 
Time-dependent expectation values of observables constructed from 
the canonical momentum are compared with thermal ones. It is found 
that a closed system, evolving from one initial condition, thermal- 
izes to high precision in the thermodynamic limit, in a time-averaged 
sense. For ensembles consisting of many members with the same en- 
ergy, we find that expectation values become stationary - and equal 
to the thermal values - in the limit of infinitely many members. Initial 
ensembles with a nonzero (noncanonical) spread in the energy den- 
sity or other conserved quantities evolve to noncanonical stationary 
ensembles. In the case of a narrow spread, asymptotic values of pri- 
mary observables are only mildly affected. In contrast, fluctuations 
and connected correlation functions will differ substantially from the 
canonical values. This raises doubts on the use of a straightforward 
expansion in terms of lPI-vertex functions to study thermalization. 



* email: aarts@thphys.uni-heidelberg.de 

* email: bonini@thphys.uni-heidelberg.de 

II email: C.Wetterich@thphys.uni-heidelberg.de 



1 



1 Introduction 



Time evolution of correlation functions in classical field theory out of thermal 
equilibrium is of interest for various reasons. On one side a classical field the- 
ory (after discretization on a lattice to regulate the Rayleigh- Jeans divergence) 
provides a relatively simple playground for testing various ideas of statistical me- 
chanics. Examples in the literature include the approach to thermal equilibrium 
[0, the relation with chaos ||, and the dynamics in thermal gradients, when the 
system is coupled to external heat baths at the boundaries ||. A more prag- 
matic reason is given by the recent emergence of classical j|, [5], |6| and semiclas- 
sical |7]] approximations to nonperturbative dynamics of low-momentum modes 
in nonequilibrium quantum field theory. Knowledge of thermalization proper- 
ties of classical field theories is important to determine the applicability of these 
approaches. 

A nonequilibrium ensemble can be defined by an initial probability distribu- 
tion with the property that the energy is distributed over the available degrees of 
freedom in some nonequilibrium fashion. One would say that an ensemble with a 
very large number of degrees of freedom has "thermalized" when after sufficiently 
long time the expectation values of observables approach the thermal ones, and 
in fact become time-independent. More precisely, we restrict this criterion to 
observables and correlation functions with support in a region of space that is 
sufficiently small as compared to the total volume. (We will always consider iso- 
lated systems, i.e. without an external heat bath.) We concentrate on equal-time 
correlation functions, and define thermal values of correlation functions by ex- 
pectation values in the canonical ensemble with the same average energy as the 
ensemble under consideration. 

It has been argued || |^, [HJ that the existence of conserved correlation func- 
tions constitutes an obstruction to thermalization. This may be illustrated by 
an ensemble consisting of many independent experiments with isolated systems 
where the initial conditions have some experimental uncertainty in the energy 
density e. As an example, this is realized in heavy-ion collisions by many inde- 
pendent scattering events. An initial nonvanishing spread Ae/e is conserved by 
the time evolution of the ensemble - it corresponds to a conserved correlation 
function. In consequence, the ensemble cannot relax asymptotically to a canon- 
ical ensemble with a temperature given by the average energy density. In the 
latter case Ae/e would have to vanish in the thermodynamic limit. In this paper 
we clarify how this issue can be related to our understanding of thermalization. 

For our investigation a natural ordering is provided by the type of ensembles 
that can be considered. Given an initial nonequilibrium probability distribution, 
we call one particular realization of this distribution a "microstate" . The sim- 
plest ensemble consists of only one microstate. We will see below that for single 
microstates fluctuations in time remain, even for asymptotically large times, and 
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we study therefore the behaviour of time-averaged observables 

(«^/' dt'Oit'). (1) 

At Jt-At 

They become stationary for sufficiently large t and At. We shall call this be- 
haviour "quasi-stationary". More advanced ensembles are formed by superim- 
posing more than one microstate: in particular, an ensemble consisting of N m 
microstates which have the same energy E and contribute with the same weight 
will be referred to as a fixed-energy ensemble. Of interest is now the behaviour 
of ensemble-averaged observables, 

1 N m 

<o>W*) = TrE<W- (2) 

Note that these expectation values can be defined locally in time, which makes it 
more suitable for the study of the time evolution of expectation values out of equi- 
librium. Expectation values in a fixed-energy ensemble with an infinite number 
of microstates are denoted as (O) e = (O) e,oo- Finally, a more general ensemble 
can be built as the superposition of states with different energy, weighted with a 
weight function f(E). 

Specific questions can now be raised. Will the nonequilibrium ensembles 
described above indeed thermalize according to the definition given earlier? Are 
there restrictions on those ensembles, e.g. on the weight function f(E)7 Is it 
important that a fixed-energy ensemble, unlike the canonical one, has no energy 
fluctuations? Further: is it possible to approximate the time evolution using a 
(truncated) expansion in, for instance, 1-particle irreducible vertex functions, at 
arbitrarily long times? 



2 Classical scalar field theory 



We will address the issues mentioned above from first principles by numerical 
simulations using a simple classical scalar field theory in 1 + 1 dimensions with 
the action 



dt / dx 



(3) 



Here L is the (one-dimensional) volume and we use periodic boundary conditions. 
The equation of motion reads d^<f> = <9^0 — m 2 — A0 3 /2, and the conserved energy 
is 



E 



J 

Jo 



dx 



(4) 



where ir(x, t) = d t <f)(x,t) is the canonical momentum. In order to solve the time 
evolution, we discretize the action in the standard way on a lattice with N spatial 
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lattice points, with spatial lattice spacing a (such that aN = L) and time step 
do, and derive the equation of motion from the discretized action.Q We then solve 
this equation numerically, starting from a nonequilibrium initial condition. 

Out of the infinite number of possible nonequilibrium initial conditions, we 
chose the following microstates: the field <j)(x, 0) is initially set to zero, and only 
a few long wavelenghts of ir(x, 0) are excited. Explicitly, 

00,0) = 0, k(x,0) =Y! Ac os{27fkx/L -^j k ), (5) 

k 

where k is integer and < ^ < 27r PJ. The prime indicates that only a very 
restricted number n c <C iV of modes is excited. We have chosen n e = 4 modes 
with momentum of the order of the mass (p = 2wk/L m m). For arbitrary 
coupling constant the initial energy is E = n e LA 2 /A. Many initial conditions 
with a certain energy can be generated by choosing different phases ipk in Eq. 
(||). We construct a fixed-energy nonequilibrium ensemble as a superposition 
of N m microstates with random ip^ and equal weight. The time evolution of 
ensemble averages can be calculated by summing over the realizations. Note that 
for N m — > oo, the fixed-energy ensemble defined in this way is translationally 
invariant: for instance the initial two-point function reads 

(tt(x, 0)7r(y, 0)) E = E'iU 2 cos[2ttA;(x - y)/L\. (6) 
k 

Microstates and ensembles with N m < oo do not have this property, of course. 

The class of initial conditions ([5]) is far from classical thermal equilibrium, 
where all modes of both <fi and 7r have nonzero expectation values of their squared 
amplitude, determined by the parameters in the action and the temperature T. 
For a free system (A = 0), the temperature is related to the total energy by 
equipartition as E = NT = LT/a. Note that this reflects the (linear) Rayleigh- 
Jeans divergence in classical thermal field theory, as a goes to zero. Corrections 
due to the nonzero coupling constant are finite in 1 + 1 dimensions. 

It is convenient to use the mass parameter m as the dimensionful scale, and 
rescale all dimensionful parameters with m, i.e. x = x'/m, t — t'/m, A = m 2 X'. If 
we also rescale the field with the dimensionless combination v = m/v3A, = v<p', 
the classical equation of motion is independent of A'. The dimensionless energy E' 
is related to unsealed energy E as E' = E/(mv 2 ) = 3\E/m 3 , and similar for the 
temperature: E'/T' = E/T. A given value of E' describes therefore several values 
of A/ m 2 if E/m is changed accordingly. Note that for fixed E/m a more strongly 
coupled theory is obtained by taking higher E'. To anticipate the Rayleigh- Jeans 
divergence in thermal equilibrium, we keep the ratio E'/N fixed, when varying 
the lattice spacing or volume. Besides E'/N, the remaining parameters that need 
to be specified are the physical volume V = mL, the lattice spacing a' = ma (or 

1 This gives the Lagrange version of the leap frog discretization. 
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equivalently N = mL/ma). For the time step we use a^/a = 0.05 (we have also 
used smaller time steps, and found no significant deviations). 

We concentrate on correlation functions that are built from the conjugate mo- 
mentum 7r(x,t), since these can easily be compared with equal-time expectation 
values in thermal equilibrium, where the only nontrivial correlation function is 

T 

(n(x,t)n(y,t)) T = -5 xy . (7) 

(X 

With respect to the space dependence, we consider both expectation values of 
global observables, i.e. observables averaged over the complete volume, and ob- 
servables that are averaged in space over a "subvolume" with size L s < L only. In 
the latter case one may think of the rest of the system as an effective heat bath, 
though it is of course still correlated with the subsystem under consideration and 
hence not a heat bath in the strict sense. 

We focus on expectation values, made from the following two building blocks: 

G?\x.,t) = ^E W - 2 (M), = ^£ ( "'V(M). (8) 

s x s x 

The label N s denotes that the average is taken over a subvolume with N s = L s /a 
points, centered around x = x s . In case that N s = N the complete volume average 
is taken, and the subscript s is omitted. The kinetic energy in a subvolume is 
given by 

K s (x s ,t) = ^aN s G?\x s ,t). (9) 

In thermal equilibrium its value is (K S ) T = N s T/2, independent of x s of course. 
One may define an effective temperature, locally in time and out of equilibrium, 
from the expectation value 

T n (t) = 2{K)(t)/N. (10) 

Note that this is an example of a "primary" quantity, i.e. it is directly given as 
an ensemble average. 

The observables we investigate are the following. In the canonical ensemble 7r 
is a gaussian (free) field, which implies that (G^)t = 3{G^)t- We want to see 
whether n becomes eventually gaussian in the nonequilibrium ensembles as well. 
We define the deviation from n being gaussian in a subvolume as 



(G^(x s ,t)) 
3(G?\x s ,t)y 



For one microstate or in a fixed-energy ensemble the conserved total energy has 
no fluctuations. The time evolution of the fluctuations in the kinetic energy of 
subsystems 

fluc(fQ = j^-y 2 . (12) 
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Figure 1: Time dependence of dev(-7r)At, averaged over mAt = 62.5, for different 
values of the energy E'/N, from E'/N = 10 (top curve, fastest relaxation) to 
E'/N = 0.5 (bottom curve, slowest relaxation). Each curve represents the average 
over 20 independent initial conditions to reduce fluctuations. Parameters are 
mL = 32, N = 256. 

is interesting for two reasons: first, although they never vanish, for a closed system 
one expects some suppression from the fact that the total energy is conserved. 
Second, one can calculate in the canonical case the size of the fluctuations exactly: 
&uc(K s ) T = 2/N s . This provides the possibility for a quantitative comparison. 
The last observable we discuss is the normalized third moment of the kinetic 
energy in a subvolume (the second moment is related to &uc(K s )), and is denoted 
as 

mom(^ s 3 ) = |g>_. (13) 

In thermal equilibrium mom(K%) T = 1 + 6/N s + 8/iV s 2 . This observable differs 
from the previous ones in that it is neither directly related to fluctuations nor 
purely local, since it involves the expectation value of {n 2 (x,t)ii 2 (y,t)ii 2 (z,t)) 
with x, y, z all in the same subvolume. These three observables are related in 
a nonlinear manner to ensemble averages, and hence examples of "secondary" 
quantities. Note also that in the canonical ensemble they are independent of the 
temperature. 



3 Single microstates 

We start with ensembles consisting of one microstate only. For each E'/N we 
generate a number of initial conditions with the same energy, and integrate the 
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Figure 2: Dependence of dev(7r)At and the relative deviation of mom(K 3 ) At from 
the thermal value on the number of lattice points N, in the quasi- stationary 
regime with mAt = 10000, for various values of mL and ma. No significant 
dependence on the subvolume size or E'/N has been found. The lines are straight- 
line fits and statistical errors are smaller than the symbol sizes. 

equation of motion. We take time averages over an interval At for each microstate 
separately. For convenience we sometimes average at the end over independent 
runs in order to reduce the fluctuations in time. 

To get some feeling for the time scales involved, we show in Fig. p] how dev(7r) a* 
relaxes towards zero as a function of time, for various values of the energy. In 
this case the time averages are taken over an interval mAt = 62.5, and since 
the size of the fluctuations in time in each microstate remains rather large, we 
show in Fig. [L] curves averaged over 20 independent initial conditions per energy 
value. As stated before, there is a correspondence between high energy and large 
coupling, and we see that the time scale involved increases for decreasing E'/N, 
as expected. 

It is clear from Fig. [I] that n tends to become gaussian: dev(7r)At — > 0. 
To make this precise, we discard the initial part where dev(7r)At is not quasi- 
stationary and calculate the time average in a long, interval mAt = 10000, 
for several values of mL and ma, subvolume sizes mL s = 4, 8, 16, and ener- 
gies E'/N = 2, . . . , 10. For each E'/N we have generated 10 initial conditions. 
We have observed no significant dependence of the long-time values on either the 
initial condition for a fixed energy, the energy E'/N in the interval indicated, or 
the subvolume size. The asymptotic average value of dev(-7r)At depends only on 
the total number of degrees of freedom N, as shown in Fig. |^. Note that changing 
mL with fixed am or vice versa has the same effect. For finite N, we see that ir 
deviates from being gaussian, on the order of a percent for iV = 100. However, in 
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Figure 3: Relative deviation of finc(K s )^ t from the thermal value, as a function 
of L s /L = N s /N, with mL s = 4, 8, 16, for various values of mL and ma. The line 
is a straight-line fit through the mL = 64, am = 0.125 points. 

the thermodynamic limit, which is iV — > oo in this classical lattice model, % be- 
comes a gaussian field. A straight-line fit shows that the thermodynamic limit is 
approached as dev(7r)At = — 1.03(2)/iV — 0.0003(2), where the numbers between 
brackets indicate the error in the last digit from the fit. Also shown in Fig. ^| is 
the relative deviation of mom(K 3 )At from the thermal value. Again we found no 
significant dependence on the subvolume size, the energy, or the particular initial 
condition. Also in this case we see a deviation from the thermal value for finite 
N, which vanishes in the thermodynamic limit as — 3.01(3)/iV — 0.0003(3). 

The relative deviation of flue (K) At from the thermal value is shown in Fig. 
|3|, again in the quasi- stationary regime with mAt = 10000, using several energies 
as before. We present the result as a function of L s /L = N s /N. For N s /N = 1, 
we see that the fluctuations are suppressed by approximately 50%, compared to 
the thermal value. A straight-line fit shows that the data are consistent with 

Zuc(K s ) At = J- (l - - 0.514(1)^ . (14) 

In other words, in the proper thermodynamic limit, N — > oo and N s /N — > 0, 
the fluctuations are given by the canonical ones. On the other hand, in the 
limit N oo with fixed N s /N, fluc(i^ s )At deviates from the canonical value, 
confirming the expectation that the heat bath must be infinitely larger than the 
subsystem under consideration. 

The results obtained so far indicate that single microstates evolving from an 
initial condition taken from the class (|j) thermalize in the thermodynamic limit, 
in a time-averaged sense. 
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Figure 4: Time dependence of dev(7r)B j jv m i n nonequilibrium ensembles with 
N m = 10 members (lower curve, shifted down by 0.2 for clarity) and N m = 150 
members (upper curve). In a larger ensemble the size of fluctuations is reduced. 
Parameters are mL = 32,N = 256, E'/N = 4. 

4 Fixed-energy ensembles 

To extend the analysis, we continue with fixed-energy ensembles, built as super- 
positions of N m > 1 microstates (|5|) with the same energy E. We present the 
effect of ensemble averaging for dev(7r) E,N m in Fig. [|for ensembles consisting of 
N m = 10 and N m = 150 members, respectively. We emphasize that there is no 
time averaging performed here, unlike the microstate case. It is seen that the 
larger the ensemble, the smaller the size of the fluctuations. To make this quanti- 
tative, we identify the size of the fluctuations with the standard deviation aE,N m 
around the mean value in the quasi- stationary regime. Explicitly, for a generic 
primary observable 

o%, Nm = {{0)l, Nm ) At - {{0)B, Nm )l t . (15) 

In Fig. | we show the asymptotic values of &E,N m for the observables T n and 
dev(7r), as functions of l/y/N m . For the quasi- stationary regime the interval 
30000 < mt < 50000 was used. We see that in an ensemble with N m — > oo mem- 
bers the standard deviations vanish, which implies that the expectation values 
become time-independent, and the ensemble becomes stationary. 

These results are consistent with the results obtained in ensembles consisting 
of one microstate only. Using definition @ for a primary observable O, we find 
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Figure 5: Size of the fluctuations <JE,N m f° r T n and dev(7r) versus l/y/N m , where 
N m denotes the number of members in a fixed-energy ensemble, ranging from 
5 to 150. The lines are straight-line fits through the eight lowest points. Same 
parameters as in Fig. [|. 

that the variance (|T5|) is given by 

4,iV ra = 4-EK°A>A*(*)-WA*(«)<Oi>A*(*)] 

m i,j 
m i 

= wD ff AM. ( 16 ) 

m i 

where cr&t,% denotes the standard deviation around the mean value in the quasi- 
stationary regime in the ith microstate. The second line follows if the individual 
microstates are statistically independent, i.e. (OiOj)At(t) = (Oi) &.t(t){Oj) &.t(t) for 
i ytz j. If we now use that the standard deviation a^t,i i n a particular realization 
is independent of the specific details of the initial condition, such that <7At,i = °"At 
for all i, we conclude that <JE,N m — &m/ 'vNmi i n agreement with the numerical 
results. Note that this implies that the initial probability distribution given by 
(||) with n e = 4 is already "rich" enough to define a fixed-energy ensemble with 
stationary asymptotic behaviour. 

5 Generic ensembles 

Perhaps surprisingly, we can use the results obtained above to show that generic 
ensembles approach noncanonical stationary distributions. Let us consider an 
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initial ensemble whose energy is not precisely fixed, but has some finite spread: 
AE 2 = ((E—(E)) 2 ) ^ 0. Note that the spread is conserved during time evolution. 
The issue of thermalization in this case is more complicated. To get some feeling 
for this, let us first consider the simple case of a superposition of fixed-energy 
ensembles weighted with a normalized gaussian function f(E): 



f(E) = A/" K exp 



1 -k{E - E) 2 /E 2 



;i7) 



with J\f K = (k/2ttE 2 ) 1 ^ 2 and where n is a large but finite constant in the thermo- 
dynamic limit. The energy spread for this ensemble is AE 2 / E 2 = 1/k. Under 
the assumption that the individual fixed-energy ensembles evolve to the corre- 
sponding microcanonical ones, the large-time asymptotic stationary values of all 
primary quantities can be obtained by folding the microcanonical ones with f(E). 
It is then straightforward to compute also secondary quantities such as connected 
functions. Using definition (|T0|) for the temperature T w , we obtainP] 



dev(7r) K « (18) 
(7r 2 (x)7r 2 (y)) K - (tt 2 (x)> k (tt 2 (y)> K « (T^/a) 2 /^ x^y. (19) 

We emphasize that secondary quantities defined for the whole ensemble are not 
weighted averages of the corresponding quantities in the individual microcani- 
cal ensembles. From Eq. (0) we see that this ensemble does not have thermal 
correlation functions, i.e. dev(7r) K ^ 0, not even in the thermodynamic limit. Fur- 
thermore, since the left-hand- side of Eq. flT9f) is nothing else than the connected 
four-point function and the right-hand-side does not vanish for \x — y\ — > oo, we 
find that the ensemble (|T7| ) does not obey the clustering property [fLl |. A viola- 
tion of clustering has a rather large impact on other observables, for instance for 
iV > N s we find that 

& ^ = w,( 1 + ^)- (20) 

Fluctuations in the kinetic energy in subvolumes will deviate substantially from 
the canonical value when k ~ N s . In the limit that N s ^> k, the size of the 
fluctuations approaches 1/ k and therefore does not vanish. 

In order to make the discussion more general, let us now consider a system 
with n c conserved intensive quantities c = {c M }, /i = 1, . . . ,n c . Supported by 
our numerical results, we assume that all microstates with given c evolve, in an 
asymptotic, time-averaged sense, to the same final state, i.e. that all details of 
the initial state other than the conserved quantities are lost. In other words, for 
a local observable O in microstate v. 

lim /' OAt) = (0)* Nc . (21) 



2 Neglecting the nonlinear T-dependence of E. 
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Here (0)* N c is the expectation value of O in the microcanonical ensemble with 
fixed c and N degrees of freedom.^ We further assume that an ensemble consisting 
of a superposition of sufficiently many states, all with the same values of c, evolves 
to the corresponding microcanonical ensemble: 

hm(0(t)) c =(0)^ c . (22) 

Notice that in this case no time averaging is required. 

We then consider ensembles that are superpositions of states with different 
values of c, weighted with a normalized weight function /aK )- As an example, 
one may consider again a gaussian distribution: 



Mc)=n 7i^r xp 



-\(c,-c,) 2 /Acl 



with Ac M ~ 1/N a . In this case Ac 2 = c 2 — c 2 ^ 0, where c = / dc f N (c)c. 
Provided that the ensemble is sufficiently "rich" , i.e. all subensembles with given 
values of c satisfy Eq. (0), the asymptotic expectation values depend only on 
the weight function /jv(c): 



(0)% = J dc f N (c)(0)* Nie . (23) 

For Ac M ~ 1/N a we extract the leading iV-dependence by writing the weight 
function fa as 

fa(c) = N an ^(s), s = N a (c - c), (24) 

where we assume that (p is in first approximation N- independent. The case a = 
corresponds to fixed Ac^/c^, independent of N. For a canonical ensemble one has 
f^ n (E/N) = Z T 1 n(E)e- E,J \ where Vt(E) is the number of states with energy 
density e = E/N and Zj< is the canonical partition function. In this case, the 
spread in the energy density is given by Ae ~ N^ 1 ^ 2 and therefore a = 1/2. In 
general, we shall assume a > 0. 

We expand (0)* Nc in a Taylor series around c such that Eq. (|23|) yields 



P) fN = E^T^T^ £~Jda<p(s)s lll ...8 



- w-^Sf^'"" 1 ' (25) 

with s^s v = J ds 99(3)5^5^. Notice that the m — 1 term in the series drops out due 
to the definition of c. In the case that is the weight function for the canonical 
ensemble, Eq. (p5|) shows that the expectation values in a fixed-energy ensemble 



5 Asymptotic values will generically be denoted with a star. 
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will differ from the canonical ones as ~ 1/JV. This is exactly the behaviour we 
found numerically (see above). Generic ensembles with a > can be compared 
with the canonical one by applying ( p5|) twice, i.e. by subtracting from (f25|) 
the same equation with = /jy"- It is then found that primary observables, 
i.e. expectation values that are directly given as ensemble averages, approach the 
canonical values in the thermodynamic limit, and the way this limit is approached 
is directly proportional to the spread in the conserved quantities For a finite 
Ac M /c M in the thermodynamic limit, the asymptotic expectation values will differ 
by a nonzero amount from both the microcanonical and the canonical ensembles, 
even in the thermodynamic limit. However, this difference is suppressed by Ac 2 . 

For secondary quantities, i.e. observables that are nonlinear functions of en- 
semble averages, the situation is more involved. For definiteness, let us go back to 
the scalar theory ([3D, which has only one conserved quantity: c\ = e.0 We restrict 
ourselves to translation-invariant ensembles. Notice also that the connected two- 
point functions for and 7r behave asymptotically as primary quantities, since 
(4>)*n = ( 7r )iv = 0- An example of a secondary quantity is the equal-time con- 
nected 4-point correlation function 

(ip(x)ip(y)ip(z)ip(w))f NiC = (ip(x)ip(y)ip(z)ip(w)) fN (26) 

- ((tp(x)ip(y)) fN (ip(z)ip(w)) fN + two perm.) , 

with ip = {ir, 0}. Here we have extended the standard definition of connected 
correlation functions to arbitrary ensembles. In particular, (7r 4 (x))j Jv c is related 
to dev(7r) in Eq. (|lT|) . Using Eq. (^5|), we find 



<- 4 (*)>Lc-(^)>W+^ 



^(^ 4 (^))W + 6 (^^ 2 ( a; ))k l 



(27) 



with s 2 /N 2a = Ae 2 . In the case that fx = (with a = 1/2), the connected 
4-point function on the left-hand-side vanishes, and the fluctuations in the energy 
density are given by s 2 = T 2 de/dT. If we use the definition of temperature in 
the fixed-energy ensemble T n /a = {tt 2 (x))* Ns and work consistently to order 1/N, 
so that T v can be identified with the temperature T in the canonical ensemble, 
Eq. (I27D yields 



(Ax))k,c=-4(l(Ax))'J —JiQ'w (28) 



and 



dev(7r) * = w^m = ~nw (29) 



The other obvious candidate, the total momentum P — ^2 x Tr(x)d x 4>(x), is not conserved 
on the lattice. We have checked that expectation values of P and P 2 approximately thermalizc, 
regardless of their initial values. 
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Figure 6: Dependence of the effective temperature T' n on the energy density 
E'/N, for several values of ma and mL. The dashed line is the free (A = 0) 
relation. 

We have measured dev(Tr) and the energy dependence of the effective temperature 
in fixed-energy ensembles independently (Figs. |2| and |], respectively), and have 
found Eq. fl29|) to be satisfied within numerical accuracy. 

As we saw in the example ([T7|) , the connected function fl26| ) does not, in 
general, obey the clustering property. Using Eq. (^), we can write 

(30) 




The limit of large spatial separation for the fixed-energy ensembles can be ex- 
tracted by comparison with the canonical ensemble, where clustering holds and 
the left-hand-side vanishes for \x — y\ — > oo. Using fx = f^ n , one finds 

hm (^ X )tf(y))* = --JL (-(tf(x))* N A +0(N^). (31) 
\x-y\-*oo ' ' N oT \de 'J 

In the microcanonical ensemble the connected 4-point function vanishes, in the 
limit of large separations, only up to terms of order 1/N. Using this result in 
Eq. (|30|) for an generic ensemble with a nonzero (noncanonical) spread in the 
energy density, we find that also in these ensembles the clustering property will 
be violated. For a = the right-hand-side of Eq. (|30"D is dominated by the first 
term on the second line and we recover the result ([TP]) . 
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It is instructive to translate these findings to momentum space. We introduce 
the notation 

(VGzi) • . .#&0) = Na6 qi+ ... +qNfi ((^i ■ ..ipn)) (32) 

where we factored out a Kronecker delta function times the volume to take trans- 
lational invariance into account. For a clustering system, all reduced connected 
functions ((...)) should be finite in the limit iV — > oo, a fixed. For the connected 
4-point function of a generic ensemble we have 

((Mnl>MYf K ,c = (tyiiM>MYf N - Na (((i>M)*f N ((i>3i>4))*f N Sq u -qi + perra.) 



N,e,C 



+ 



d 2 



2N 2a de 2 



(tyiiMMYNie, 



c 



s 2 a 



N 2, 



Q-l 



^{^MYn^^MYn^Ax -?2 + Verm. 



(33) 



When momenta come in pairs, the reduced connected 4-point functions for the 
microcanonical and canonical ensembles differ even in the thermodynamic limit, 
and, remarkably, for ensembles with a < 1/2 these functions actually diverge. A 
similar calculation for the connected 6-point function gives, to leading order, 



s 2 a 2 



N 2 ° 



d 



(34) 



d 



((M2)Yn^{{^)Yn^((^MYnA 



' de 



de 



"92^93, 



-94 



+ perm. 



A comparison with the canonical ensemble shows that the reduced connected 
6-point function in the microcanonical ensemble diverges in the thermodynamic 
limit ~ N. For higher order connected functions to be finite for all momentum 
combinations, has to agree with f™ n to higher and higher order in 1/N. Ap- 
proximation methods based on expansions in connected or 1-particle irreducible 
correlation functions may break down since higher-point functions could strongly 
affect the expansion. Without taking care of the failure of clustering properties 
they may not be suitable to describe thermalization phenomena. This might ex- 
plain some qualitatively different results obtained in ||, 10 , where the same scalar 
theory was studied using truncations of an exact flow equation for the time evo- 
lution of equal-time correlation functions |J. It also puts doubts on the validity 
of approximations which lead to thermalization for arbitrary initial ensembles. 



6 Conclusion 

Let us summarize our conclusions for the time evolution of classical ensembles 
with nonequilibrium initial conditions, based on numerical results for the <p A the- 
ory in 1 + 1 dimensions on the lattice. The investigated statistical ensembles 
approach in the large-time limit stationary ensembles that depend only on the 
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probability distribution for the energy density (and the other intensive conserved 
quantities). For ensembles that consist of a small number of microstates, the 
asymptotic ensemble is approached only in a time-averaged sense. For a very 
large number of degrees of freedom (thermodynamic limit) we find that ensembles 
with a nonzero energy density spread do not approach the canonical ensemble: 
their large-time primary observables differ from the thermal values by a finite 
amount proportional to the spread. Ensembles with vanishing spread, on the 
other hand, do thermalize, up to terms suppressed by the number of degrees of 
freedom. These corrections can be computed, and our numerical results were 
found to be in excellent agreement with the analytical results. 

We also notice that generic statistical ensembles do not cluster for large time. 
Secondary quantities (like connected or 1PI functions) receive corrections that in 
coordinate space do not vanish in the limit of large spatial separation; as a result, 
corrections to quantities defined over large subvolumes can be very large. In mo- 
mentum space, violation of clustering implies that for specific momentum values 
secondary quantities diverge in the thermodynamic limit. Clustering requires a 
vanishing initial spread in all conserved quantities with extreme precision. We ar- 
gue therefore that care should be exercised when applying approximate methods 
based on expansions in connected or 1-particle irreducible correlation functions 
to the study of thermalization: such methods might in fact break down in the 
large-time regime. 
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